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A combined experimental and theoretical study of very- 
low-energy electron diffraction at the (0001) surface of 2H[- 
NbSe2 is presented. Electron transmission spectra have been 
measured for energies up to 50 eV above the Fermi level with 
k" varying along the FK line of the Brillouin zone. Ab initio 
calculations of the spectra have been performed with the ex- 
tended linear augmented plane wave k-p method. The exper- 
imental spectra are interpreted in terms of three-dimensional 
k"-resolved one-electron band structure. Special attention is 
paid to the quasi-particle lifetimes: by comparing the broad- 
ening of the spectral structures in the experimental and cal- 
culated spectra the energy dependence of the optical poten- 
tial —iVi is determined. A sharp increase of Vi at 20 eV is 
detected, which is associated with a plasmon peak in the 
Im[— 1/e] function. Furthermore, the electron energy loss 
spectrum and the reflectivity of NbSeg are calculated ab ini- 
tio and compared with optical experiments. The obtained 
information on the dispersions and lifetimes of the unoccu- 
pied states is important for photoemission studies of the 3D 
band structure of the valence band. 



I. INTRODUCTION 

Very-Low-Energy Electron Diffraction (VLEED) spec- 
troscopy with incident electron energies below ~ 40 eV 
has recently been established as the experimental method 
giving a direct access to dispersions E{k) and lifetimes 
of quasi-particle states above the vacuum level (see 
Refs. |^-|4 and references therein). The VLEED spec- 
tral structures reveal the critical points in the surface- 
perpendicular dispersions E{k^) of the states that couple 
to the incident electron beam to form the LEED state. 
The fundamental advantage of the VLEED spectroscopy 
over the conventional unoccupied band structure meth- 
ods, such as inverse photoemission or X-ray absorption 
spectroscopy, is that it does not involve initial states and 
in particular cases, for example at high-symmetry lines 
parallel to the surface, is capable of restoring 3D Bloch 
vectors k. __ 

Within a quasi-particle picturcfi'El the time-reversed 
LEED states serve as the final states of the one-step 
photoemission theory. Independent information on their 
dispersions and lifetimes provided by the VLEED ex- 



periment can be used in the photoemission experiment 
to resolve the valence rtftnd structure in the surface- 
perpendicular directionJj'EI which is often blurred by the 
complicated non-free-electron-like dispersion and strong 
self-energy effects in the final states. In particular, the 
broadening of the VLEED spectral structures provides 
information on the lifetimes of the final state quasi- 
particles, which can be used to estimate the fc"*" broad- 
ening of the Bloch constituents of the LEED state and 
thereby to judge on the intrinsic uncertainty of the band 
mapping. 

Within a simplified one-particle approach, the LEED 
problem reduces to solving the Schrodinger equation for 
a semi-infinite crystal given the energy E and the initial 
conditions of the incident electron, i.e., the asymptotics 
of the wave function in the vacuum. In the plane parallel 
to the crystal surface the LEED wave function obeys the 
Bloch theorem and is characterized by the 2D Bloch vec- 
tor kll. In the vacuum, far from the crystal surface, it is a 
superposition of plane waves: the plane wave propagating 
towards the crystal defines incident current, and the to- 
tal current carried by the LEED state is the transmitted 
current. The ratio of the two currents is the transmission 
coefficient T{E). In the bulk the LEED wave function is 
a superposition of propagating (real fc"*") and evanescent 
(complex k^) Bloch waves. 

This approach neglects inelastic scattering, so that in 
the case of a non-zero transmission the LEED state neces- 
sarily includes a propagating Bloch wave. In what follows 
the propagating constituents are referred to as conduct- 
ing states. Gross features of the T{E) spectrum depend 
upon the band structure of the bulk crystal: the ener- 
gies at which the band ceases transmitting the current 
(e.g. when the group velocity vanishes) and other criti- 
cal points are reflected in experimental spectra. For ex- 
ample, T{E) drops abruptly to zero when an energy gap 
in the k" projected band structure is encountered. How- 
ever, owing to inelastic processes, one never observes zero 
transmission in the experiment, and instead of a sharp 
drop predicted by the simple theory one observes a rather 
smooth decrease of the transmission as a sign of the crit- 
ical point. J-. 

In 1937 Slateiu showed that the broadening of the spec- 
tral structures and the absence of energy gaps can be re- 
produced by adding an imaginary term, the optical po- 
tential —iVi, to the potential in the crystal half-space. 
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Slater's idea is to associate the effect of the optical poten- 
tial with a spatial damping of the wave functions rather 
than with a decay in time. Then the Bloch vector ac- 
quires an imaginary part of the order of Vi/Tiv^ (u^ is 
the group velocity), or, in terms of the mean free path 
d = 1/Imfc, it is = fiv^/d. The optical potential 
is, thus, understood as the inverse lifetime of the quasi- 
particle. Alternatively, one can think of an electron with 
a lifetime r = Ti/Vi and mean freejjath l/Imfc that moves 
in a real non-absorbing potential.EJ In the present work we 
make an attempt to extract quantitative information on 
the energy dependence of Vi from the measured T{E) 
spectra by comparing the broadening of experimental 
and calculated spectral structures. 

The unoccupied band structures of layered quasi- 
2D materials, giWitL, as graphiteEl or transition metal 
dichalcogenides,0'oEj are particularly interesting: ow- 
ing to the interlayer potential barrier, a nearly-free- 
electron model fails to provide even a qualitative pic- 
ture of the unoccupied states. The necessity to describe 
such states within an all-electron approach for the crys- 
tal potential of general shape and to understand their 
relation to the diffraction process has been a driving 
force of the development of the Bloch- wave-based ah ini- 
tio methodology.llirtil In this work, we use the extended, 
linear augmented plane wave (ELAPW) k • p method,t3 
whose distinctive feature is that it reduces the inverse 
band structure problem - finding k-^ given k" and E - 
to a matrix eigenvalue problem. 

Of the layered materials, the band structure of 211- 
NbSe2 deserves special attention. Its Fermi surface com- 
prises two cylindrical-sheets and a tiny 3D electron pocket 
near the F pointJlB'tJ all having different superconduct- 



ing properties.E3 Mapping of this pocket in the pho- 
toemission experiment is complicated by that its ex- 
tension in 1^1— is comparable with the final state 
broadcningJijnJ and it relies critically on the knowledge 
of not only the final state dispersion but also lifetime. 

In this paper, we present experimental electron trans- 
mission spectra of NbSe2 and compare them to our ah 
initio calculations. The experimental technique is de- 
scribed in Sec. ^ A brief account of our self-consistent 
band structure calculation is given in Sec. [II and the 
computational method for LEED states is presented in 
Sec. [V. The experimentally observed broadening of 



spectral structures is reproduced by including an imag- 
inarypart of the potential VJ into the Hamiltonian. In 
Sec. |y| we describe a procedure to obtain quantitative 
information on the quasi-particle lifetimes from the anal- 
ysis of the shape of the experimental spectra. The ob- 
served energy dependence of the optical potential is in- 
terpreted in terms of electron energy loss function. In 
Sec. we present the unoccupied band structure of 
NbSe2 in the FK direction and discuss possible limita- 
tions of the common description of the inelastic scattering 
by an optical potential. Optical properties are discussed 
in Sec. Ivil. 



II. EXPERIMENTAL 

Essential details of ourjJZLEED experimental tech- 
nique are given elsewhere. BEil Briefiy, we used a stan- 
dard 4-grid LEED optics operating in the retarding field 
mode. In this mode the electrons are accelerated in 
the gun to the energies required to form a well-focussed 
beam (normally from 100 to 300 eV) and then deceler- 
ated in a retarding field between the gun and the sam- 
ple. This mode allows us to achieve the lowest pri- 
mary energies without any significant degradation in fo- 
cussing. However, angle dependent measurements are 
complicated by that the retarding field distorts the off- 
normal electron trajectories. To determine K", we pa- 
rameterized its dependence on energy E and the sample 
rotation angle a by a biquadratic function K^\E^a) = 

2 

^ AimE^a™ with the coefficients Aim fixed by fitting 

i,m=0 

to the experimental points with well-defined if", namely 
to the angle dependent target current onsets, in which 
K^\ = ^J2m {E — e4>) /h, and to the characteristic diffrac- 
tion patterns when a diffracted beam hits the electron 
gun exit, in which K" equals half the surface reciprocal 
lattice vector GII/2 (see Ref. ^ for details). 

The present experiment was carried out on an im- 
proved setup: Screening from the stray magnetic fields 
using /i-metal shielding reduced the displacement of the 
electron beam on the sample to less than 0.5 mm over 
the whole energy range. The position of the diffracted 
beams in the area obscured by the electron gun was con- 
trolled using a miniature fiuopjscent screen covered by 
ZnO:Zn low-energy phosphor,LJ which was mounted at 
the gun exit. Careful adjustment of the incidence angles 
was achieved using a custom made sample holder with 
3 angular degrees of freedom. This is particularly im- 
portant for the fiaky crystals of layered materials, whose 
surface after gluing to the sample holder can appear a few 
degrees off. The tilt angle adjustment was coupled to the 
linear shaft motion of our standard manipulator, and the 
azimuthal angle was adjusted using an additional fork 
mounted on a rotary wobble stick. This design guaran- 
teed the incidence angle adjustment accuracy better than 
0.25° in all angular degrees of freedom. A symmetric 
design of the sample holder prevented distortion of the 
electron trajectories in the retarding field. The energy 
spread of the primary beam was ~0.25 eV HWHM. The 
size of the beam spot on the sample can be estimated as 
less than 0.5 mm through the whole experimental energy 
range. 

As a practical matter, it should be noted that during 
the angle-dependent measurements the electron spot can 
slightly displace along the surface due to the retarding 
field and mechanical drifts in the manipulator. This dis- 
placement can be a problem for the layered materials, 
whose crystals closer to the periphery typically incorpo- 
rate minor misoriented crystallites. The sample position 
should therefore be optimized to avoid the misoriented 
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crystallites in the whole experimental energy and angle 
range. This can be controlled most easily by inspecting 
the diffraction patterns on the screen. In this respect the 
VLEED experimental setups based on the LEED electron 
optics have an advantage over the inverse photoemission 
based ones. 

The experimental system provided two measurement 
modes: 1) measurements of the elastic reflectivity into 
individual diffracted beams, often referred to as the 
VLEED technique, using the LEED screen and a CCD 
camera, and 2) measurements in the target current, com- 
monly referred to as the Target (or total, absorbed) Cur- 
rent Spectroscopy (TCS) technique - see, e.g., Ref. |2^. 
The VLEED mode gives most detailed information about 
the diffraction process, but its principal limitation is that 
the reflected intensities in the area obscured by the elec- 
tron gun can not be measured (e.g. the specular beam 
near the normal incidence). The TCS mode is free of 
this problem and also benefits from greater experimental 
simplicity, but its limitation is that the target current 
I{E) gives the total reflectivity which includes the elas- 
tic reflectivity R{E) integrated over all diffracted beams 
and, in addition, the inelastic reflectivity Rinci{E) corre- 
sponding to the secondary electrons that leave the crys- 
tal. The inelastic contribution gives, however, only a 
rather featureless background, so the structures in the 
I{E) curves reflect essentially the elastic electron trans- 
mission T{E) = I — R{E), which contain the band struc- 
ture information (the derivatives dl/dE and dT/dE are 
practically equivalent). 

The term VLEED will further be used in reference to 
the dominant physical mechanism forming the spectral 
structures in the individual beams as well as in the total 
reflectivity, rather than to the experimental technique. 
In this context we will refer to the target current spectra 
also as the VLEED spectra. 

Atomically clean (0001) surfaces of NbSe2 were ob- 
tained by standard cleavage in the vacuum chamber at a 
base pressure of 7 x 10~^° mbar. Compared to other lay- 
ered materials, the lifetime of the NbSe2 surface is small: 
in a day after cleavage the spectral contrast significantly 
decreased and the background in the LEED patterns in- 
creased, indicating strong surface contamination. 



III. BAND STRUCTURE 

The calGHiations are performed with the ELAPW k • p 
method.tBESH The self-consistent potential was con- 
structed within the local density approximation (LDA) 
of the density functional theory with the full-potential 
augmented Fourier components technique described in 
Ref. |2^. The basis set included 487 energy independent 
APWs (energy cutoff 10.2 Ry), and the extension of the 
radial basis set contributed another 200 basis functions. 
The extension was introduced following the prescriptions 
of Ref. p^. The Brillouin zone (BZ) integrations were 



performed by the tetrahedron method with a mesh of 
550 k points that divides the irreducible BZ into 2187 
tetrahedra. 

All the occupied states down to the Nb 4p semi-core 
band (at -31 eV) are treated as valence states. The 
density-of-states (DOS) function is in good agreepient 
with the FLAPW calculation of Corcoran et alx3 and 
wiiii X-ray photoemission measurements of Wertheim et 
aVE3 The five peaks due to the strong Se Ap - Nb Ad 
hybridization have distinct counterparts in the experi- 
mental spectrum, see Fig. ^ The peak in the unoccu- 
pied DOS at 2.7 eV arises from localized states of almost 
pure Nb d character. Thus, one can expect strong dipole- 
allowed transitions from the Nb Ap ban d to produce an 
absorption peak at fiLO ^ 34 eV (see Sec. VII). We do not 
observe any localized states above E — Ep — 5 eV; the 
complicated fine structure of the DOS curve for higher 
energies reflects a non-free-electron-like behavior of the 
unoccupied states. 



IV. CALCULATED VLEED SPECTRA 

Our calculations are based j-an the Bloch waves ap- 
proach to the LEED problemj23 in which the crystal is 
treated as a semi-infinite system and the LEED wave 
function in the crystal half-space is sought as a linear 
combination of propagating and evanescent solutions of 
the Schrodinger equation for a bulk crystal potential. 
Our implementation of the theory in the case of a sin- 
gular all-eiectron crystal potential has been described 
elsewhere.EJ Now we briefly sketch the computational 
procedure. 



A. Inverse Band Structure Problem 

In the ELAPW-k • p method the complex band struc- 
ture can be obtained by an analytical continuation of the 
Schrodinger equation to the complex k spaccEj In ap- 
plication to semi-infinite systems, we use the ELAPW- 
k • p method to solve the inverse band- structure problem: 
given two real Cartesian components of the Bloch vector 
k" = {kx, ky) and the energy E, wc find the values of k^ 
that satisfy the Schrodinger equation 



iJ*(i;,kll +zfc,-|:;r) =i;*(£;,kll +-z,k^,]r) 



(1) 



for the Bloch vector with a complex z-component k^ (z 
is a unity vector in the z-direction) . 

In the k • p method the wave function of a Bloch vector 
(k'l -t-zfe-*-) is a product of a trial function with the Bloch 
vector (k'l -|-zfcj]-) and the phase factor exp[i(fc"'^ — fc,^)z]. 
The trial function is a linear combination of the ba- 
sis functions ^j(ko;r) constructed at the reference point 
ko = k'l -I- 2,k^ with a real k^: 



^'(£;,kll+zfc^;r) = e 



[EM) ^,(ko;r) 
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inverse problem can be written as a matrix equation 



H + 25iP^ + {5i -E)d 
S; 



Cn = 0, 



(2) 



In terms of the coefficient vectors C„ = {C,''" } the Here rll = {ra:,ry) and ci' are surface reciprocal vectors. 

To construct the Laue representation we use a rapidly 
converging 3D Fourier decomposition of the wave func- 
tion, which differs from the true all-electron wave func- 
tion only in a close vicinity of the nucleus (see the goug- 
ing technique described in Refs. ^ and ^9|). Each sur- 
face Fourier component fsn{z) is smoothly continued into 
the vacuum half-space by a linear combination of two 
plane waves, one of which is the outgoing vacuum solu- 
tion exp [z/csz] with | kH -|- ci' p + fc^ = E, and the other 
one is a decaying (unless ci' ~ 0) plane wave exp [i^sz] 
with a purely imaginary z-component Qs- In the present 
calculation it has been chosen k^ — = AE for all s 

with AE = 3 Ryd. For GjJ = one of the waves repre- 
sents the incident electron wave, and the other one the 
reflected wave, i.e., qo = —k^. 

Thereby the functions are defined in the whole 
space and they are smoothly continuous everywhere. The 
LEED function $ is sought as a linear combination of the 
functions that minimizes the value 



with H being the Hamiltonian, O the overlap matrix, 
and the z-projection of the momentum matrix. By 
introducing the vector i5„ — — {2P^ + ^n^) we re- 
duce Eq. (^ to a linear eigenvalue problem of twice the 
dimension: 



H-EO 

1 2P-L 



= s, 



I 

-6 



Dn 



(3) 



Thus, the solutions of the inverse band-structure problem 
k:^ are obtained as the eigenvalues S:^ of this generalized 
non-Hermitian problem. 

Because of the rather large unit cell of NbSe2 (six 
atoms per unit cell) the solution of Eq. (^), espe- 
cially for complex matrices, results in prohibitively time- 
consuming calculations .Ea It has turned out that the CPU 
time can be greatly reduced by orthogonalizing the orig- 
inal basis set, i.e., by seeking the trial function in terms 
of the bulk Hamiltonian eigenfunctions '!/)j(ko; r) at kg 
rather than in terms of the basis functions ^j(ko;r). 
Eliminating the overlap matrix alone makes the comput- 
ing of Eq. (H) three times faster. The price for that is 
the necessity of recalculating the momentum matrix , 
which scales quadratically with the number of both ba- 
sis functions and eigenfunctions. To further accelerate 
the procedure we cut the eigenvalue spectrum and re- 
tain only tpj with energies lower than 5 Ryd above the 
spectral interval we are interested in. A suitable cut- 
off energy is difficult to determine a priori, however, the 
quality of the results can be verified by comparing the 
resulting inverse band structure k'^{E) with the usual 
band structure E{k^). 



B. Constructing the LEED Function 

In the bulk half-space the LEED function <& is ex- 
panded in terms of the solutions = '^{E, k" -f z/c^; r) 
of Eq. (||). The quality of the wapj functions is known 
to deteriorate with growing 

IC I -c^ With the present com- 
putational setup reliable solutions were found to be re- 
stricted to the interval ImJ^j; < 3 A"-'^. To construct the 
LEED function, we take all the solutions in this inter- 
val, which for NbSe2 is typically ten (at the vacuum level 
Evuc) to twenty Bloch waves (at 40 eV above i?vac)- Then 
a Laue representation of the functions is constructed 

*„(rll,z)= ^ /,„(z)exp[z(kli+Gli)rll]. (4) 

s=0 



{H -E)<^\\ 



\{H ~ E)^{y)\^ dv, 



(5) 



under the constraint that the incident current is equal to 
unity. The resulting function ^(r) is the closest to a solu- 
tion of the Schrodinger equation in the surface region. By 
construction the trial function satisfies the Schrodinger 
equation both in the bulk half-space and far from the sur- 
face in the vacuum, so that only a finite slab in the vicin- 
ity of the surface contributes to the integral (||). This 
slab may include the part of the crystal where the poten- 
tial strongly deviates from the bulk potential, but in this 
calculation we adopt a simplified step-like shape of the 
surface barrier, so the potential changes abruptly in the 
middle of the van-der-Waals gap from its bulk distribu- 
tion to the constant vacuum value. The generalization of 
the method to the case of an arbitrary potential in the 
interface region has been described in Ref. |30| . 

The minimization procedure 5\\{H — E)^\\ = leads 
to a system of linear equations. If the auxiliary tails 
exp [ i(7s-z ] accidentally cancel then the result is equivalent 
to the simple matching of wave functions. It is, however, 
difficult to predict whether the matching is possible and 
to separate the matching error from the errors introduced 
at thfi stage of obtaining the LEED function constituents 
.113 If the exact matching is impossible the variational 
method yields a solution that is smoothly continuous by 
construction - this function is thought to be the best ap- 
proximation for the LEED function we seek. The accu- 
racy of the ultimate result can be estimated a posteriori 
by checking the current conservation in the LEED func- 
tion, i.e., by comparing the current transmitted into the 
crystal by the superposition of the propagating Bloch 
waves to the current in the vacuum carried by the 
superposition of plane waves. 
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C. Normal Incidence Spectrum 

The k'l — spectra along with the underlying real 
band structure are presented in Fig. ^. The Vi = calcu- 
lation yields partial currents carried by the Bloch states 
(shown by whiskers). It reveals the conducting states 
responsible for the transmission of the current, thereby 
offering an interpretation of the VLEED spectrum in 
terms of conducting fragments of the band structure: the 
abrupt changes in T{E) all reflect critical points in the 
conducting bands. 

The deviation of the k^{E) points by the inverse 
ELAPW-k- p from the E{k-^) lines by the direct method 
with the full basis is seen to be negligible. Together with 
the reasonably small current non-conservation this justi- 
fies the basis set reduction we have undertaken to make 
the calculations feasible. 

A special feature of the normal incidence spectrum is 
that the conducting bands do not overlap, i.e., there is 
only one dominant propagating constituent in the LEED 
function. In this respegt the picture is similar to that 
of the Slater ID model.El The inclusion of the imaginary 
term lifts the principal difference between the propagat- 
ing and evanescent waves (whose properties to couple 
to the incoming wave may be very similar) and thereby 
levels out the intensity variations (see the dashed curve 
in Fig. 1^). In particular, already at a moderate value of 
Vi = 0.5 eV the narrow gaps in the T{E) spectrum almost 
completely disappear, so one can hardly expect them to 
be observed in the experiment. Wider gaps remain well 
visible, although the reflected intensities are strongly re- 
duced and the structures broadened with respect to those 
predicted by a Hermitian Hamiltonian. In the next sec- 
tion we will use these properties of the Slater theory to 
extract quantitative information about the damping of 
quasi-particles from the shape of experimental spectra. 

V. DETERMINING THE ENERGY 
DEPENDENCE OF Vi 

We shall now compare our measured target current 
spectra I{E) with the calculated transmission coefficient 
T{E). First of all, we need to bring the I{E) spectrum 
to the same absolute units as T{E). This is not triv- 
ial because of the presence of the unknown background 
caused by the secondary electrons contributing to the 
reflectivity. In view of the close similarity of the fine 
structure of the measured I{E) and theoretical T*'^(i?) 
spectra over a wide energy region, see Fig. 0, we have 
chosen to determine the experimental T{E) curve by 
fitting I{E) to T'^^{E) with the linear transformation 
T <"=?(£') = aI{E) +b + cE. Here the function b + cE 
represents a linearly varying background. 

Strong intensity variations in the T{E) spectrum are 
much more important for the baad structure information 
than the shape of the maxima,Q which are very broad 



even if the damping is neglected (see Fig. ||). Thus, for 
a conclusive comparison, the extrema of the dT{E)/dE 
function must be analyzed. 

Owing to computational instabilities, theoretical T{E) 
curves contain numerical noise, which should be excluded 
before the derivative is calculated. To remove the jagged- 
ness we average the T{E) function within an interval of 
width Ai? by fitting it to a parabola: for a given E a 
parabola is constructed that approximates the raw T{E) 
spectrum over the interval [E — lS.E/2] E + AE/2] and its 
value at E is taken to be the intensity of the smoothed 
T{E) spectrum. The width AE of the interval is 0.6 eV 
at the vacuum level and it grows linearly up to 2.4 eV at 
50 eV above -Evac- Thereby slow changes of intensity are 
not affected and only the noise is removed. 

We extract the energy dependence of the optical po- 
tential Vi{E) from the experimental spectra by compar- 
ing the sharpness of the dT/dE extrema (maxima and 
minima) in the experimental and theoretical spectra. To 
quantitatively characterize the sharpness, we associate 
the function dT/dE in the interval between its two zeros 
El and E2 with a parabola that passes through zero at 
El and E2 and embraces the same area S as the dT/dE 
function, S = T{E2) - T{Ei). Then the curvature of the 
parabola 65*/ {E2 — Ei)^ is taken to express the sharpness 
of the structure. This parameter changes very strongly 
over the spectrum, so in Fig. ^ we show the cube root of 
the curvature \^S/ {E2 — E1). We have performed a series 
of calculations for the normal incidence with an energy 
independent Vi ranging from 0.5 to 3 eV and considered 
the sharpness of each structure as a function of optical 
potential. Three examples are shown in the right panel 
of Fig. ^. The value of sharpness from the left panel is 
taken to yield the value of the optical potential in the 
right panel. One immediately notices a sharp increase 
of the optical potential between 16 and 21 eV: the two 
minima change their shape very similarly with increasing 
the Vi (circles and diamonds in right panel of Fig. ^) and 
the minimum at 21 eV in the experiment is significantly 
broader than the one at 16 eV (see left panel of Fig. 

The values of Vi determined in this way are neces- 
sarily overestimated. First, in addition to the inelas- 
tic processes in the electronic system, there exist broad- 
ening mechanisms dependent upon experimental condi- 
tions, such as surface roughness and finite energy res- 
olution. Also the angular spread of the incident beam 
contributes to the broadening owing to the kH dispersion 
of spectral structures. Secondly, our theoretical spec- 
tra do not take into account inelastic scattering in the 
surface barrier region, e.g., scattering on the surface de- 
fects. This effect results in an additional reduction of 
the spectral structures but, unlike the inelastic scatter- 
ing in the bulk, it hardly affects their energy broadening. 
These mechanisms should, however, be less significant for 
the layered materials because of the rather small electron 
density at the surface and a small concentration of defects 
on the cleaved surface. On the other hand, our absolute 
values of VJ are in accord with the recent ab initio results 
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on the j-quasi-particle lifetimes for noble and transition 
metals:Ed For energies around E — Ep = 5 eV the calcu- 
lated lifetimes do not exceed a few fs {Vi ~ 0.2 0.7 eV). 

The energy dependence Vi{E) derived from the experi- 
ment is presented in Fig. ^. To get an idea of how reliable 
the obtained Vi (E) dependence is we compare the normal 
incidence data (circles) with the data for k" — jTK (> 
(The off-normal incidence spectra are shown in Fig. 
i|rK| = 0.3 A). The absolute values of Vi agree very 
well and both spectra suggest a sharp increase at around 
20 eV. The non-steady behavior of the k" = ^FK derived 
values around 30 eV seems to have no physical meaning; 
it reflects the disagreement between theory and experi- 
ment in shape of the T{E) maximum at 25 eV, which in 
section VI we explain by the deviation of the calculated 
band structure from the experimentally observed one. 

An optical potential is associated with the imaginary 
part of the electron self-energy E, whose energy depen- 
dence is expected to reflect singularities of the energy 
loss function — Jm[l/e(q, w)]. In particular, in the GW 
approximation,Ej Im E is given by an integral, whose in- 
tegrand contains the inverse dielectric function (DF). We 
are not in a position to calculate the lifetime of the LEED 
constituents ab initio. A rough idea of the average effect 
of plasmon excitations on the energy dependence of the 
electron lifetime can be obtained from the DF for q = 

E 

by plotting the integral / — Im[l/e(q — 0,uj)] duo. We 



have calculated the DF in the random phase approxima- 
tion, see Sec. VII, In Fig. || we compare the integral func- 
tion to the experimental dependence Vi{E). The sharp 
increase of the optical potential agrees well with the plas- 
mon location at 21 eV (see also Fig. |l^ in Sec. VII ), 
which is confirmed by the trpismission energy loss mea- 
surements by Bell and Liang.LJ Our theoretical DF shows 
also a second plasmon peak at 35 eV, which may be re- 
sponsible for the step-like growth of the optical potential 
at higher energies. 



VI. TK DIRECTION 

Experimental and theoretical T{E) curves for kll rang- 
ing from zero to FK are presented in Fig. ^. The exper- 
imental and theoretical spectral structures show a good 
agreement in both energy dispersion and relative ampli- 
tudes. Owing to the increase of the incident beam spread 
in k'l , the low-energy part of the experimental spectra at 
large k" looks smeared compared to the theory. Note 
that because of the vector k" changing with the inci- 
dent electron energy the spectra cannot be compared in 
a one-to-one manner. To facilitate the comparison of the 
dispersion of the spectral structures we present a gray- 
scale r(kll,E) plot. Fig. ^. A very good agreement be- 
tween the measured and calculated k' ' dispersion of the 
T{E) intensity proves that the theoretical approach we 
have adopted is adequate. The results have turned out 



rather sensitive to details of the potential distribution in 
the bulk half-space; we stress that the advantage of the 
ELAPW-k • p method of accurately taking into account 
the strong non-mufhn-tin effects makes the method an 
indispensable tool for studying the electron diffraction 
on layered materials. 

We now analyze the information on the bulk band 
structure contained in our experimental data (and de- 
liberately blurred by the optical potential in the theory) . 
The critical points in the fc"*" dispersion of the conduct- 
ing bands (e.g. the band edges) manifest themselves as 
the extrema in the dT/dE curve. Of course, reflected is 
the optical-potential-damped (completely complex) band 
structure, in which the band dispersions are smoothed, 
and the critical points arC|-shifted by some tenth of eV 
from their Vi = positionsJd The deviation of the surface 
barrier from the step-like form may also cause shifts of 
dT/dE extrema, but for the layered material this effect 
seems to be insignificant. 

To stress the k" dispersion of the spectral structures 
we construct a gray-scale plot. Fig. ^ in which the shad- 
ing shows the energy area covered by the conducting 
bands (the regions between dT/dE maxima and min- 
ima); k'l projected band gaps then appear as white areas. 
Again one observes excellent overall agreement between 
the experiment and theory. The plot, however, reveals 
also differences in finer details, which reflects fundamen- 
tal limitations of the theory employing a one-electron 
band structure and optical potential. The deviations 
are necessarily present due to our neglect of the real 
part of the self-energy E(£',k) - it is replaced by an 
exchange-correlation potential in the local density ap- 
proximation. For example, the white area ranging from 
kli =Q,E-Ef = 27eVtokll = 0.7 A"!, 31 eV in the ex- 
periment is shifted towards higher energies by about 1 eV 
with respect to its calculated counterpart. In particular, 
this causes the higher width of the measured T{E) max- 
ima between 25 and 29 eV than in the theory (see the 
discussion in Sec. 0). The difference is especially well vis- 
ible between 0.2 and 0.4 A""'^. The discrepancies are seen 
to significantly increase above ^40 eV: an overall shift of 
the experimental bands to higher energies is observed. 

The deviations are seen to increase at larger k" . From 
the band-structure point of view the main difference be- 
tween the normal and off- normal incidence spectra is that 
in the latter case the conducting bands intersect; in other 
words the LEED state contains more than one propagat- 
ing Bloch constituent. An example for k" = FK is shown 
in Fig. ^, lower panel. The theoretical T[E) curve for 
negligible absorption (solid line) has a different charac- 
ter from the normal incidence spectrum (Fig. H): there 
are no energy gaps in the band structure along the KH 
line, and the transmission never drops to zero and never 
reaches unity - it almost never exceeds 0.8. Thus, not 
only minima but also maxima are affected by the op- 
tical potential. Also the positions of spectral minima 
strongly depend upon the function Vi{E). For example, 
with increasing Vi, the minimum at 19.6 eV moves away 
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from its measured location towards lower energies and 
in the spectrum calculated with the energy dependent 
it appears at 17.9 eV. Such a high sensitivity of results 
to a phenomenological parameter may be an important 
limitation of the theory. Although the overall shape of 
the k" = FK spectrum is well reproduced by the the- 
ory, the band structure information may be incorrectly 
transferred to the T{E) curve. 

Presumably, the weak point is the presence of several 
propagating Bloch waves in the band-structure decom- 
position of the LEED function. In connecting the optical 
potential to the electron attenuation we basically rely 
on the relation Imfc^ ~ Vi/v^. The parameter Vi is, 
thus, easy to operate when there is only one wave, but 
its meaning is not that transparent when there are sev- 
eral waves with different velocities (which, in addition, 
are connected by the matching conditions). This point 
of view is supported by Fig. |l^ in which we compare the 
measured k" dispersion of the T(E) minima with theoret- 
ical results for the energy dependent Vi and for a constant 
moderately small Vi of 0.5 eV. The latter results are ex- 
pected to provide undistorted information on the band 
structure. One can see that whenever there is a disagree- 
ment between small Vi and large Vi results, for small A:'! 
the experimentally determined function Vi{E) brings the 
positions of the minima closer to the experiment (areas 
A and B). On the contrary, for large /c" the Vi = 0.5 eV 
results are often closer to the experiment (areas C and 
D). 

Such deviations can be corrected by introducing a po- 
sition dependent (in general non-local) 14, which would 
allow to control the behavior of the complicated LEED 
state. In this case different Bloch constituents would be 
damped differently depending on their spatial distribu- 
tions. Such an approach within the multiple-scattering 
formalism of VLEED has been developed by Bartos et 
al. in Ref. El 



VII. OPTICAL PROPERTIES 

Optical properties of the crystal provide independent 
information on the unoccupied band structure. Within, 
the one-particle theory - random phase approximation^ 
without local field effects - the calculations employ the 
same u noccu pied bulk states as the LEED calculations 
of Sec. fVAl. 



The imaginary part of the dielectric function e2(w) was 
calculated as an integral over the BZ with the tetrahe- 
dron method (for computational parameters see Sec. III). 
The real part ei(w) was determined out of e2(w) by 
the Kramers-Kronig integration with the energy cutoff 
^(^max = 70 eV [The values of £2(1^) above 38 eV are un- 
derestimated because of the finite number of bands taken 
into account - from Nb 4p (-31 eV) up to 38 eV above 
the Fermi level.] 

In Fig. ^ we compare our ah initio results with the 



reflectivity measurements by LiangESI and, transmission 
energy loss spectrum by Bell and Liang.E3 The experi- 
mental data are in excellent agreement with the calcula- 
tions. The main plasmon peak at 21 eV corresponds to 
the termination of transitions from the Ad valence band 
states. Our calculations predict also a rather strong peak 
at 36 eV due to the transitions from the semi-core Nb Ap 
states to the unoccupied Nb Ad peak dX E — Ep = 3 eV 
(see Fig. |l|). The high-energy plasmon is important to 
understand the energy dependence of the electron life- 
time. 



VIII. CONCLUSIONS 

Angle-dependent VLEED measurements in the target 
current mode have been performed on the layered 211- 
NbSe2. Tight control over the electron beam and adjust- 
ment of the incidence angle in three angular degrees of 
freedom enabled highly accurate measurements for the 
incidence wave vector k" scanning the FK line of the 
Brillouin zone. Based on the self-consistent band struc- 
ture, ab initio VLEED spectra have been calculated with 
the ELAPW-k • p method, the inelastic scattering being 
modelled via the optical potential. 

A novel methodological aspect is the algorithm to ex- 
tract the information on the electron lifetime from the 
VLEED experiment. The obtained results are consistent 
with available experimental and theoretical information 
on the dielectric function of NbSe2 . 

Owing to the complicated non-free-electron-like unoc- 
cupied band structure of NbSe2, the VLEED spectra 
show rich structure over the energy region up to 50 eV 
above the Fermi level. A good agreement between ex- 
periment and theory in the energy location of the TCS 
structures over the whole interval makes it possible to 
interpret the observed shape of the spectra in terms of 
the energy dependent optical potential. 

The comparison of the experimental and theoretical 
results focusses on the surface-projected band structure 
represented by the k" dispersion of the extrema in the 
energy derivative of the transmission coefficient dT/dE. 
A good agreement between the experiment and theory 
proves that our theoretical approach is capable of accu- 
rate description of the unoccupied bands of NbSe2 . Mi- 
nor discrepancies have been observed, which are traced 
back to the imperfectness of our one-electron approach as 
well as to intrinsic shortcomings of the phenomenological 
treatment of inelastic scattering. 
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FIG. 1. Density of states for NbSe2: total DOS (solid 
line) and /-projected Nb d (dot-dashed line) and Se p (dotted 
line) partial DOS per atom. (Note the energy scale break at 
7 eV.) Dashed curvp-phows the X-ray photoemission spectrum 
by Wertheim et alx 



FIG. 2. Calculated unoccupied band structure of NbSe2 
in the FA direction. The lines in the left panel show the 
band structure E{k^) for Vi — obtained with the direct 
ELAPW-k ■ p method (without the basis reduction). The 
conducting bands are marked by whiskers, whose upper ends 
show the band structure k'^{E) obtained with the inverse 
method. The length of the whisker is proportional to the 
contribution of the Bloch wave to the target current. Trans- 
mission coefficient T{E) (calculated in the vacuum half-space) 
is shown in the right panel by the solid line. The shaded area 
shows its deviation from the current in the bulk (the sum of 
the group velocities times the weights with which the con- 
stituents enter the LEED function). The dashed line is the 
T{E) spectrum with the optical potential Vi — 0.5 eV added 
to the Hamiltonian H in Eq. (H). 
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FIG. 3. Comparison of the normal incidence target cur- 
rent I{E) and transmission coefficient T{E) spectra of NbSe2. 
Tlie measured I{E) curve is shown by the dot-dashed line. 
Theoretical T{E) and dT/dE spectra are shown by dashed 
lines in the lower and upper panels respectively; they are 
calculated with an energy dependent optical potential [the 
Vi{E) function is shown in Fig. |5|. To arrive at the experi- 
mental T{E), the original I{E) function is fitted to the theo- 
retical T{E) by a linear transformation that involves scaling 
and subtracting the background. The resulting experimental 
T{E) and dT/dE curves are shown by solid lines. 

FIG. 4. Left panel: sharpness of the extrema in the ex- 
perimental (x) and theoretical (o) dT/dE spectra (see upper 
panel of Fig. H) . Negative values are ascribed to minima and 
positive to maxima. Right panel: dependence of the sharp- 
ness upon the optical potential for two minima (at 16 and 
21 eV) and a maximum (at 30 eV) in the dT/dE spectrum. 

FIG. 5. Dependence of the optical potential Vi on the 
incident electron energy extracted from the normal (o) and 
an off-normal (x) incidence spectrum. The function Vi(E) 
used in the calculations is shown by the dashed line. The 

E 

solid line is the function J — Im[l/e(q = 0,aj)] dio. 





FIG. 6. Experimental (left) and theoretical (right) angu- 
lar dependence of the T{E) spectra with the vector k" scan- 
ning the FK line of the BZ. In the experiment the sample 
rotation angle a increases in steps of 0.5°. The onsets of the 
spectra are at the energy etp + {h^ /2m)K^^ at which the elec- 
trons start penetrating into the solid. Owing to the retarding 
field, the incidence angle differs from a and depends on en- 
ergy in such a way that k" varies along the spectrum almost 
linearly. In the left panel the k" values are indicated that cor- 
respond to the onset of the spectrum and to _E — Ep ~ 50 eV. 
Note the scale change at 18.5 eV (marked by a dashed line). 



FIG. 8. Experimental (left) and theoretical (right) k" 
dispersion of the dT/dE spectra. [Original T{E) curves are 
shown in Fig. ^. The shading fills the regions between the 
maxima and minima in the dT(E) / dE spectra, i.e., the area 
where the second derivative d^T/dE^ is negative. Physically, 
the shaded area shows the surface-projected dispersion i5(k") 
of the conducting bands. The gray scale shows the amplitude 
of d^T/dE^ (in a logarithmic scale), which characterizes the 
amplitude and sharpness of the extrema. The shown energy 
region begins 0.5 eV above the transmission onset. 

FIG. 9. Dependence of the electron transmission spectrum 
for k" = FK on the optical potential. Experimental spectrum 
is shown by the dashed line. In the lower panel conducting 
fragments of the band structure in the H-K-H interval are 
shown (see also the caption of Fig. S) 



FIG. 10. k" dispersion of T{E) minima: experiment 
(dots); theory with the energy dependent Vi (circles) and with 
Vi = 0.5 eV (lines). 

FIG. 11. Optical properties of NbSe2 for two light polar- 
izations: E ± c (solid lines) and || c (dashed lines). The 
experimental reflectivity spectruma for E ± c is shown in the 
upper panel by the dot-dashed line. The triangles in the lower 
panel show the energy locatji«n and intensity of the low-energy 
structures of the measuredE3 EELS spectrum. The fragment 
of the Lorentzian fiimction (dotted line) reproduces the exper- 
imentally observed^ energy location, intensity, and width of 
the main plasmon. 



FIG. 7. Experimental (left) and theoretical (right) k" dis- 
tribution of the T{E) spectral intensity shown in a linear gray 
scale. White areas correspond to maximal T{E). [Original 
T{E) curves are shown in Fig. ^. The energy region shown 
begins 0.5 eV above the transmission onset. Note the scale 
change at 18.5 eV (marked by a white line). 
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